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Abstract: We discuss the dependence of the results of renormahzed perturbation the- 
ory for dilute polymer solutions on the choice of the uncritical manifold where the 
perturbation series are evaluated. Special emphasis is given to the influence of polydis- 
persity corrections on the results of one and two loop calculations. For monodisperse 
solutions we establish that after a Borel resummation the dependence on the choice of 
the uncritical manifold decreases when higher orders of the expansion are taken into 
account. 
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1. Introduction 

In now more than thirty years since de Gennes' observation pp that the correlations of a 
single long polymer chain can be mapped on the critical correlations of a m-component 
ferromagnetic spin model in the limit m — 0, the application of dilatation symmetry to 
dilute polymer solutions has led to a rich and mature body of knowledge, which is able 
to explain the universal features found in experimental data. These developments have 
been reviewed in the book of des Cloizeaux and Jannink [2j and with special emphasis 
on excluded volume effects in the monography of Schafer 

Since exact results for realistic models are available only in space dimension = 2, 
quantitative calculations for the physically interesting case d = 3 heavily rely on the use 
of perturbation series expansions. Unfortunately, the expansion parameter diverges in 
the critical limit of polymer length n —>■ oo. This can be remedied by mapping the result 
of perturbation theory on a uncritical manifold of the parameter space with the help of 
the renormalization group (RG). The RG maps the physical variables chain length n, 
monomer size / and excluded volume strength (3^. onto a set of renormalized variables 
{'nR, Iri u), where low order perturbation theory for a given scaling function can be safely 
evaluated yielding reasonable results. The determination of the uncritical manifold, by 
choosing the renormalized segment size Ir introduces some numerical parameters into 
the theory, which we will fix in the sequel by fitting our results for universal ratios to 
values measured independently in simulations or experiments. Once those parameters 
are fixed, the theory has to prove its quantitative accuracy by making predictions for 
additional measurable quantities. The necessity of such a fit procedure is clearly due to 
the error which we introduce by calculating the perturbation expansion only to low order. 
With this approximation we break the strict scale invariance of the full renormalized 
scaling function under the RG map and introduce a dependence of the results on the 
choice of the renormalized length scale Ir. We expect that the influence of the choice of 
the uncritical manifold will vanish gradually when more and more orders of perturbation 
theory are taken into account. One main goal of the present paper is to check this 
expectation for those few observables, namely the mean square end-end-distance Rl, the 
radius of gyration and the second virial coefficient A2, where the perturbation series 
have been pushed to six, four and three loop order respectively. Since most of the more 
complex observables of interest have been evaluated only at zero or one loop level, we 
spent some effort to describe the optimal choice of the renormalized manifold at the one 
loop level as discussed in [3] . We also discuss how one can handle polydispersity effects 
within the theory, a important topic when we try to explain experimental measurements, 
which never work with purely monodisperse samples. 

Our paper is organized as follows: In section 2 we define our polymer model and the 
observables of interest. We also review some basic notations of polydispersity. In section 
3 we review the RG map and discuss the choice of the renormalized manifold at the 
one loop level. In section 4 we present the results of higher order perturbation theory 
and check their dependence on the choice of the uncritical manifold before and after a 
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suitable Borel resummation. Section 5 gives a conclusion of our findings. A collection 
of perturbative results is presented in the appendix. 

2. Polymer Model and Observable Quantities 

We represent a polymer chain of n segments in a simple spring and bead model by n + 1 
beads, linearly connected to their neighbors by elastic springs with mean distance I and 
interacting with each other via a repulsive (5-pseudo-potential of strength (5^. This leads 
to the hamiltonian 

. n+l 

^{^^■> = E(^^- - ^^-0' + i^^iy^ E ^(^^- - ^r)^ (2-1) 

J=l j,3' 

where Vj is the position vector of the bead number j and d denotes the space dimension. 
The partition function 

V[r]e--^^^^ with V[r]^\{j^^ (2.2) 

for dimension d > 2 can be calculated only in a perturbation expansion which orders 
in the parameter z = /^e'T-^"'^^^- Similar considerations hold for averages of observables, 
defined as 

<0>^^ J V[r] O e-^{''^>. (2.3) 

In the sequel we focus our interest on the mean square end-end-distance 

Rl =< (fn - fof > (2.4) 
and the radius of gyration 

j=0 

where Rem — Z]j=o '^i center of mass vector of the molecule. The second virial 

coefficient can be read off from the virial expansion of either the osmotic pressure U 

U = c, + ^Afcl + 0{4), (2.6) 

where Cp is the polymer concentration, or of the forward scattering intensity 

1 

cr\q = 0, c„ TV) = — + -Ac + 0{c% (2.7) 

where c = CpN is the monomer concentration and A^^ is the weight averaged chain 
length defined as below. The two definitions of and coincide for monodisperse 
systems but differ for a general chain length distribution P[n). Both quantities can be 
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obtained from the second virial coefficient A2{ni,n2) for two chains of lengths ni and 
n2 according to ^3] by averaging over the chain length distribution P{n)\ 

^^=Y. P{ni)P{n2)A2{ni,n2) (2.8) 
^2 = E ^^K)^(^2)^K,n2), (2.9) 

ni,n2 

where = ^„^-P(^) is the average chain length. Two other standard chain length 
averages that show up in the literature are the weight average A^^ and the 2;-average 
defined as 

N^--=^ry,P{n)n^ = Np2 and := ^y^P{ny = . 

(2.10) 

A^^ and can be expressed as indicated above in terms of the average chain length 
and the second and third moments p2 and of the reduced chain length distribution 
p{y) defined by 

P{n) = ^P(^)- (2.11) 
3. Renormalization 

In order to map our perturbative results to the uncritical manifold we introduce 
renormalized variables according to 

I =XIr (3.1) 
n = X^'^riaZniu) (3.2) 
/3e = X'uZM, (3.3) 

where e = 4 — d and the scahng parameter AG [0, 1] measures the degree of dilatation. 
For finite polymer concentration Cp and finite momentum q we define renormalized 
quantities via CpR = Cpl\ and qR = qlu. The Z-factors Z„ = and Zy_ = ^ have been 
calculated for 0^ field theory in the minimal subtraction scheme to five loop order (see 
|3] and refs. in there) and are given in the appendix for e = 1. The dependence of the 
Z-factors on a change of the scaling parameter A can be obtained via integration from 
the flow equations 

= W{un) (3.4) 
X^^\n& = 2-^ (3.5) 

X-^HZ^) =viu), (3.6) 

where the Wilson function W{u), except for a linear dependence on e, and the exponent 
functions r] and u depend only on the renormalized coupling u. Since the perturbation 
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series for W, rj, and u are only asymptotic, they have to be resummed to yield reliable 
results. We follow the work of Schloms and Dohm [S] who resummed the flow functions 
at the upper critical dimension dc = 4: and then evaluated the flow equations and 
renormalized scaling functions directly in d = 3 dimensions, without further expansion 
in e. Besides the Gaussian fixed point at u = the Wilson function W{u) has a nontrivial 
fixed point at u* = 0.364, which is related to the excluded volume limit n —>■ oo, j3e > 0. 
The correlation length exponent u, which governs the power law ~ N'^'^ for R"^ and 
i?g in the excluded volume limit, takes the fixed point value z/(w) = 0.588. To measure 
the distance from the excluded volume fixed point, we introduce the parameter f = ^■ 
We can now use our perturbative results for the observables Re,Rg and to form 
universal ratios, i.e. combinations depending only on / and global characteristics of the 
system like space dimension or polydispersity. Such quantities reduce to pure numbers 
at the fixed points of the RG. Just like the critical exponents, they are universal in the 
sense that they are independent of the microstructure of the underlying model. From 
our three observables we can form two independent ratios 

QRl ^ f d 



g 



where the prefactor of the interpenetration ratio if)^ , which roughly measures the volume 
that a chain excludes for other chains, has purely historical reasons. For polydisperse 
systems, where Af and differ, another ratio can be introduced by using A^. 

Choice of the uncritical manifold in one loop approximation 

As a general recipe, the renormalized length scale Ir should be chosen smaller than 
the smallest macroscopic length scale important for the observable of interest. For 
example, it does not make sense to choose Ir > Rg, because then the whole coil would 
be smaller than one effective segment of size Ir. In lowest order perturbation theory we 
find Rl = —I'^N, which after renormalization gives 

Rl = llN.R. (3.8) 

Thus the choice N^r ~ 1 fixes Ir as Ir ~ Rg. Note that with this procedure we 
introduce some polydispersity dependence into the choice of the uncritical manifold. 
We will discuss this later in the section. Now finite polymer concentration Cp and finite 
momentum q both introduce additional characteristic length scales into the system j3] 
and depending on their values and the observable of interest it can be necessary to choose 
~ or //J ^ i in the appropriate limits. Following [3] we choose the relation 

qt> N,R Co 

to interpolate smoothly between the above limits. The constants ^q, Uq and Cq have been 
introduced in ()3.9p . since our qualitative arguments fix Ir only up to a constant and the 
dependence of our results on the parameters qQ,no and cq displays the approximation 
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we made by truncating the perturbation expansion at low order. In the sequel we will 
concentrate on the dilute low momentum limit Cp — > 0, g 0, where ()3.9|) reduces to 
NzR = no- For readers interested in the determination of and cq we refer to [Sj. 
Following we use the fixed point value of the interpenetration ratio = ^^^^^u*) to 
fix no, since it depends on N^r already in the zero loop approximation. 
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Figure 1. Zero and one loop results for the universal ratio ^P* plotted versus uq: The 
upper curves display the result for a monodisperse ensemble and the lower ones for a 
exponential ensemble. The dashed line shows the one loop result for an exponential 
ensemble evaluated with the choice Nji = Uq. 

For a monodisperse ensemble {p2 = 1) and an exponential ensemble {p2 = 2) the 
zero and one loop approximations are displayed in figure^together with the most precise 
value = 0.247 obtained from monodisperse computer simulations jH] which is in fair 
agreement with the best experimental value = 0.245 ±0.005 [7]. While the zero loop 
result reproduces the simulation result for a value uq = 1.85 the one loop result allows the 
choices uq = 0.53 and tiq = 2.95 to fit the simulation data. In j^] the value no = 0.53 was 
chosen mainly because of the unreasonably large polydispersity corrections that occur in 

when the relation Npt = no is used for the determination of the uncritical manifold. 
We can resolve this problem by absorbing the polydispersity dependence partially in the 
choice NzR = no (corresponding to Ir ^ Rg) of the uncritical manifold. For the one loop 
result of the exponential ensemble both choices of the uncritical manifold are included in 
figure — reffigl. One finds that our choice Nzr = no allows us to choose no = 2.95 with 
a value for ^* = 0.112 consistent with the zero loop result 0.11 for rio = 1.85 and close 
to the value 0.12 found with = no = 0.53 in 0. Unfortunately presently neither 
experimental nor simulation data on the polydispersity dependence of \1/ are available. 
We prefer the choice no = 2.95 since it enhances the numerical precision of several one 
loop results for intra chain properties. As an example we consider two universal ratios 



6 



aR and 5 which have been studied in the hterature on di-block copolymers. We divide 
the chain into two blocks of relative length xi = ^ and X2 = which may have 
different chemical compostition. This setup allows for two different values uu and U22 
of the intra-block excluded volume repulsion and for a third value U12 of the inter-block 
excluded volume repulsion. The ratios (Tr and 5 are defined as 

T>2 ,02 and 5= ' ' ' , (3.10) 

^elH + ^e2H ^^l^-iK^glH + ^g2H) 

where the subscript H denotes zero inter-block coupling U12 = 0. At the symmetrical 
fixed point U12 = uu = U22 = u* the ratios and S can be evaluated explicitly. Using 
the asymptotic power law i? ~ n'^ we find jH] 

For other values of the excluded volume strength we can evaluate an and 6 only 
perturbatively and the result again depends on the choice of the renormalized manifold. 
The radii of gyration of both blocks are additional length scales which have to be 
considered in the choice of the renormalized manifold. In the definition of and S 
one observes that in the limit Xi —>■ the contribution of the smaller block Xi can be 
neglected compared to the contribution of the larger block. Thus we can safely use 
NizR + N2zR = Uq for the determination of the uncritical manifold, as in the homo- 
polymer case. 




Figure 2. One loop and exact results for the universal ratios ct^j and S* for equal 
block size xi — X2 = ^ plotted versus no 



Figure |21 displays the no dependence of the one loop results for cr and S for the 
symmetric case xi = X2 = ^, together with the exact values of eq. p.lip . One finds 
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that the choice no = 2.95 reproduces the exact values within a accuracy of 2%. In figure 
El we plotted the exact and one loop results for the fixed point values and 6* as a 
function of the relative block size xi. The perturbative results deviate from the exact 
ones (using u = 0.588 from high order perturbation theory) by less than 2% over the 
whole interval of block compost it ions. Similar results have been found for a one loop 
calculation of the persistence length Lp jH]. 

The limitations of the one loop approximation can be judged from the result for the Uq 
dependence of the universal ratio R"^^^ displayed in figure El For hq = 2.95 the value 
-^g/e ~ 0.983 is still closer to the value 1 for noninteracting chains than to the value 0.96 
obtained in high precision simulations of self-avoiding chains at the excluded volume 
fixed point [TUj . 

1.30 T 




Figure 3. One loop (dashed) and exact (full lines) results for the universal ratios a'^ 
and 6* for uq — 2.95 plotted versus xi 

4. Higher orders of perturbation theory 

For only a few polymer observables, namely Rg, Re and A2 the loop expansion has been 
pushed beyond the two loop level. Nickel and coworkers calculated Re up to the order 
of six loops jTH 1121 to ^^^^ loops |T3j and A2 to two loops for a monodisperse en- 

semble. Their results enable us to investigate the uq dependence of the universal ration 
-^g/e t*-* ^'^^^ loop level. In order to further study the behavior of we calculated 
the three loop contribution to the second virial coefficient A2{n,n) ^1], which allows 
the evaluation of at the three loop level for a monodisperse ensemble. In addition 
we calculated the two loop contribution to R^ and \I' for arbitrary polydispersity. 
First we want to test wether a consistent choice of no is possible at the two loop level. In- 
specting figurelHwe find that the two loop result for ^P* reproduces the value = 0.247 
for a value of uq = 4.78. Furthermore one observes that the uq dependence of is 
rather weak for a range of values no G [0.5,4] but at a value of ~ 0.18 — 0.2 well 



8 



below the desired result. Using the condition N^r = Uq leads to \E'* = 0.085 for the ex- 
ponential ensemble, somewhat below the zero and one loop results. The choice Nr = uq 
(dashed curve) instead does not allow for a consistent fit with reasonable polydispersity 
corrections. With the choice uq = 4.78 we can read off the prediction R"^*^ = 0.963 from 
the two loop result in figure El a value that compares already well with the Monte-Carlo 
result 0.96 and with the two loop epsilon expansion result 0.959 fTB] . 
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Figure 4. Two loop results for ^* plotted versus no- The upper full curve corresponds 
to a monodisperse ensemble. The lower full and dashed curves correspond to a 
exponential ensemble evaluated with N = no and Nji = no respectively. 
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Figure 5. One to four loop results for the universal ratio R^*/„ plotted versus uq. 
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Figure 6. One to three loop results for the universal ratio plotted versus hq. 
Beyond the two loop level the three and four loop results for Rj*, displayed in figure 

9/ ^ 

EJshow the formation of a plateau at the value 0.963 for hq e [2, 6], nursing the hope that 
the choice of no becomes less important with increasing order of perturbation theory. 
On the other hand our three loop result for in figure El exhibits a pronounced no 
dependence that seems to contradict our expectation. We can trace back this behavior 
to the asymptotic nature of the renormalized perturbation series, which have to be 
resummed in order to extract sensible information beyond the two loop level. 



Resummation 

It is well known that the perturbation expansions in quantum field theory usually 
are only asymptotic series with zero radius of convergence jTH]. This gives rise to an 
exponential growth of the expansion coefficients. The leading behavior of the coefficients 
in high order perturbation theory can be obtained from a semiclassical calculation as 

mniHi 

Pk = k\i-afk'c + , (4.1) 

where in d = 4 the coefficients are a = | for our definition of the coupling and 
b = 2 + M for a correlation function involving M polymer chains. With the knowledge 
of the asymptotic behavior we perform a standard Borel resummation procedure of our 
perturbation series for R^, Rg and A2 as described in jTHllin]- The coefficients bo > 
and a involved in the resummation procedure have been tuned to 60 = 6 and a = 1 
in order to give optimal convergence of the approximations. Note that despite the fact 
that we evaluate the perturbation series directly in d = 3 dimensions we were forced to 
use the d = 4 result a = | in order to obtain good convergence of the resummed series. 
This may be traced back to the fact that for the renormalization we used Z-factors 
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which have been defined via minimal subtraction of e poles at dimension d = 4(e = 0). 
These Z-factors also are given as asymptotic series and the stronger exponential growth 
of their coefficients seems to dominate the asymptotics of the renormalized series at 
d = 3. A similar procedure was used in the calculation of universal quantities of the 
0{m) symmetric 0^ model j^. 
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Figure 7. Resummed one to three loop results for the universal ratio plotted 
versus uq. 




Figure 8. Resummed one to four loop results for the universal ratio R^*/^ plotted 
versus no- 

Figures [7| and |H1 display the resummed results for \1'* and R'^*^, obtained via ()3.7p 
from the resummed functions for Rl,Rg and Af. For both ratios the variation as a 
function of Hq is greatly reduced - the range of uq in figures d and |H1 is extended by a 
factor 5 compared to figures El and IHl The most prominent effects can be seen in the 
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two and three loop results for The plateau region, where the two loop result was 
insensitive to uq is shifted close to the value expected from simulations. The variation of 
the three loop result is greatly reduced, being now fairly insensible to no in the interval 
Uq e [2, 30] for a value around 0.247. The effects of the resummation on i?^*^ are less 
dramatic. Mainly the value of the plateau already present in figure El is shifted from 
0.963 to 0.96 and thus is now in full accord with the simulation results. 



5. Conclusions 



We studied the dependence of several universal ratios on the choice of the uncritical 
manifold, where the evaluation of renormalized perturbation theory gives sensible 
results. We found that the inclusion of the polydispersity dependence of the zero loop 
result for the radius of gyration into the choice of the uncritical manifold allows for the 
choice no = 2.95 leading to reasonable polydispersity corrections and to improved one 
loop estimates for several universal ratios. Furthermore we find this procedure necessary 
in order to obtain a consistent choice of rio at the two loop level. It would be interesting 
to check our predictions on the polydispersity dependence of by comparison with 
precise experimental or simulation data. Beyond the two loop approximation we found 
that a resummation of the asymptotic series is mandatory to establish the increasing 
insensitivity of the results on the choice of no, which was expected on theoretical grounds. 
The resummed series for and Rf/^ display extended plateau regions where the no 
dependence is weak, at values \E'* = 0.247 and Rfje = 0.96 in full agreement with 
experimental measurements and Monte-Carlo simulations. 

Appendix A. Perturbative results 

Bare monodisperse results 

The perturbation series for R^, R^ and A2{n, n) evaluated directly in li = 3 dimensions 



are 



Rl = 6i?2(i + - y)/ + 6.2968797/ - 25.0572507/ 

+ 116.134785/ - 594.71663/), (A.l) 
R, = Rla + + - + 6.564897/ - 26.70629/), 

(A.2) 

A, = (47r)3/2i?3^(l - 32(7 -4v^) ^ ^ 13.92783/ - 80.30/), (A.3) 

X5 



where R?. = and z = P^n^^'^. 
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Z-factors 

The Z-factors as obtained from (j)'^ field theory [4^ evaluated for d = 3 are 

Zn = l-u-^u^- 1.2708333m3 - 5.299419m^ + 40.504065M^ (A.4) 
1 43 

Zu = -(1 + 4m + —u'^ + 43.639293m^ + 7.439240u^), (A.5) 

Zi = l--u- —u^ - 0.91666667^^ - 3.266246159M^ (A.6) 
2 16 

Copolymer quantities 

The Z-factors for the copolymer case at one loop order evaluated for d = 3 are 120] 

4""'Hwaa') = ^ (1 + + t^a'a') + 2Uaa' + 0{u')) (A.7) 

Z'l^\uaa) =1-Uaa + 0{U^). (A.8) 

The direct evaluation of renormalized perturbation theory in = 3 dimensions leads to 

1 /2 

Rla = lln.n{l + Uaa {^^/J^ " l) + ^ {3U.'J^ + US.'J^ 
^ 384Kt^ + 640kI^ + 176kI^ - 32k,^ + 13 \ \ 

(1 + /€ai?)^/2 J J 

for the radius of gyration of block a, where KaR = denotes the renormalized ratio of 
the length of both blocks and to 

(l + ^^^^f^, ((67 + 196.,.).;/? - (105 + 315.,.); 
+ io5(i''^'^^^)3 ((67 + 196.2fi)4f ^ - (105 + 315.2i?)) 

+ 105(1 + k,rY V^^^ + "^^^ - ^^^^ + '^^^^ + '^^^ ) 



- 196(1 + k^rY'\kir + .^/j') + 268(.i. + kI^) + 402.?^)^ 

(A.IO) 

for the radius of gyration of the whole chain. For the renormalized mean squared end- 
end distances we find 

1^,, C^r,^/^ 1^^2?iaa I/2A ^ 3/2 8.^^ + 4/t,. - 1 



(A.ll) 

/ 2^1/2 _ ^ 2^1/2 _ 

i?2 = 6/^(^1. + 712.) 1 + ^^11 + ^22^^ 

V 1 + 1 + ^ 



+ l^2R 



13 



+ 



2u 



12 , 



IriiR + n2R 



,1/2 



1 I 3/2 
(1 + KanY'^ 



(A. 12) 



Renormalized polydisperse results 

For a general chain length distribution p{y) we obtain the following two loop result for 
the second virial coefficients A2 and , where m — corresponds to A2 and m = 1 to 
The polydispersity correction vanishes for a monodisperse ensemble. 

- ^fNl (1 + «7 (2 - VN~nf^ (7 - 4v^) ) + ^z* 
V 15^ ^ ^ 4 V 15 12 



1024^2 327r , 125 

\ In 9 J 



15 



In 2 + 
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arctan 
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^ ^ dy y"'+^ p{y) + 



■Pm+1 Jo 



Pm+1 Jo 



(A.13) 

dy p{y) 
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^00 roc 

dyi / dy2 yTvTiyi + 2/2)^ Piyi)piy2) -7 + 4^2 



Pm+1 Jo 

+ [2(^-1) -?|v^^ , 
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dy y"*+t 
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The radius of gyration for a polydisperse solution can be obtained from the small 
momentum behavior of the density correlation function This leads to the following 
average of the radius of gyration Rg{n) of an isolated chain 



RI[P]= r dy^-^y'R\yN) 
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(A.15) 
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(A.16) 

1 

where the coefficients ak,Zu{k) and Zn{k) are taken from ()A.2|) . ()A.5|) and ()A.6|1 
respectively. Again the polydispersity correction Cg vanishes for a monodisperse 
ensemble. 
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